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ABSTRACT 


The  computational  method  developed  in  this  thesis 
permits  the  calculation  o-f  the  aerodynamic  performance  of  a 
wing  o-f  arbitrary  plan-form.  Both  basic  and  additional  1  1  f  t 
are  analyzed.  This  treatise  is  restricted  to  thin  wings  in 
steady,  inviscid,  i ncomprsssi bl e  -flow.  The  method  uses  a 
grid  system  o-f  control  points  over  the  winq  semi -span.  The 
circulation  over  the  wing  is  considered  variable  with 

discrete  values  at  the  speci-fied  grid  points.  Finite 

/  : 

difference  equations  are  utilized;  to  determine  these 
discrete  values.  Control  point  i ndetermi nac i es  are 
evaluated  analytically.  Matrix  inversion  is  required  for 

solution  by  the  method  presented.  Details  of  the  matrix 
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technique  are  developed  in  Ref.  4.  A  brief  summary  of  the 
pn  ncipal  computational  relations  is  included.  No 

numerical  results  are  yet  available  but  are  expected  during 
the  next  phase  of  this  research.  ' 
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BACKGROUND 


The  overall  research  of  which  this  thesis  forms  a  part 
has  a  three-fold  purpose.  Firstly,  it  is  intended  to 
instruct  the  student  about  the  fundamental  equations  that 
govern  the  aerodynamics  of  wings  in  incompressible  flow. 
Secondly,  it  teaches  the  student  how  to  convert  the 
fundamental  equations  into  a  form  suitable  for  numerical 
calculation  on  the  digital  computer,  and  to  carry  out  the 
complex  programming  required  for  this  purpose.  Thirdly, 
once  the  first  two  ob j ect i ves  have  been  ac comp  1 1  shed ,  the 
computer  program  itself  provides  a  useful  pedagogical  and 
design  tool  for  illustrating  the  effects  of  various  wing 
design  paramenters  on  the  final  aerodynamic  performance  of 
a  wi ng . 

Lcdr .  J .  L .  Par ks  made  a  creditable  start  toward  the 
first  abjective  in  his  thesis  CRef.  13.  However,  in  the 
limited  time  available,  he  was  unable  to  male  significant 
progress  toward  the  second  or  third  objectives. 

It  was  evident  early  in  this  present  analysis  that 
Parks'  computational  technique  had  some  flaws.  The 
expectation  was  that  these  flaws  would  be  readily 
discovered  and  soon  corrected  and  that  objectives  two  and 
three  would  be  rather  quickly  and  easily  completed.  Tnat 
has  not  been  the  case1 


R 


As  the  investigation  actually  developed,  a  whole  series 
ot  obstacles  wss  encountered  and  each  was  in  turn 
eventually  overcome.  The  progress  was  neither  quid  nor 
easy.  In  every  case,  it  was  -found  that  the  seeming 
obstacle  was  really  attributable  to  some  conceptual  error. 
Once  the  error  was  -found,  the  obstacle  disappeared  and  the 
1 n vest  1  gator s  gained  a  deeper  insight.  Some  of  the  major 
obstacles  that  were  encountered  and  eventually  overcome  are 
discussed  in  the  following  paragraphs. 

Confusion  about  the  proper  type  of  computational  grid 
to  use,  whether  staggered  or  unstaggered,  and  the  grid  size 
necessary  for  accuracy  was  a  major  factor.  An  8  by  8 
unstaggered  grid  over  the  semi -span  is  now  used,  for  good 
r eason . 

Another  point  of  confusion  was  how  to  deal  numerically 
with  the  singularity  that  occurs  in  the  governing  integral 
equation.  One  option  is  to  avoid  the  problem  by  staggering 
the  grid  of  field  points  with  respect  to  the  grid  of 
control  points.  The  other  option,  which  was  chosen  -for 
this  investigation,  resolves  the  indeterminacy  by  rigorous 
analysis  and  as  a  consequence  allows  the  employment  o+  a 
simple  unstaggered  mesh. 

A  third  point  of  con+usion  concerned  the  evaluation  of 
the  partial  derivatives  of  the  circulation  function. 
Should  analytical  or  finite  difference  methods  be  used  to 
represent  these  derivatives^ 


n 


It  was  found  that  analytical 


G  i  f  f  sr  ent  1  at  1  on  is  incorrect  and  that  -finite  differences 
must  be  used. 

Still  another  confusing  ooint  concerned  the  validity  of 
representing  the  circulation  function  by  a  Fourier  series. 
While  this  procedure  is  widely  advocated  in  the  technical 
literature,  it  was  found  that  a  much  simpler  and  clearer 
formulation  can  be  obtained  otherwise. 

Boundary  conditions  also  caused  some  confusion. 
Special  conditions  apply  at  the  leading  and  trailing  edges 
and  at  the  wing  tip  and  midspan.  These  quite  cample;; 
boundary  conditions  have  now  been  fully  and  rigorously 
analysed  and  incorporated  into  the  formulation  of  the 
or obi sn. 

It  is  evident  from  the  foregoing  discussion  that  tms 
research  has  amounted  to  a  major  education  in  basic 
aerodynamics  and  in  numerical  methods.  In  these  respects 
it  has  been  a  richly  rewarding  experience. 

However,  in  view  of  the  foregoing  obstacles  and 
problems,  the  time  schedule  has  of  course  been  great! / 
delayed  from  that  which  was  initially  anticipated.  Thus, 
there  are  no  final  numerical  results  at  this  particular 
stage  of  the  investigation.  There  is  also  no  more  time 
available  to  this  investigator.  Hence,  this  final  aspect 
wi ! 1  have  to  be  completed  by  some  subsequent  investigator. 
Nevertheless,  what  the  present  effort  nas  produces,  to  aio 
anv  subsequent  worker,  is  a  very  sophisticated  and  retired 


ana  cne 


Again,  reversing  tne  limits  of  integration 


3.  Iqsorciic  sign  the  inteprand,  tn?  lntepral  scjua^  1  c-n  for 
tne  inauceo  velocity  due  to  the  trailing  vortic: ty  Decodes 


D.  TOTAL  VORT I Cl  TV  EFFECT 

The  total  induced  velocity  at  a  particular  control 
point  is  the  summation  of  the  effect  produced  oy  the 
vorticity  over  the  wing  and  the  effect  produced  ov  the 
trailing  vorticity.  This  means  that  for  each  control 
point,  the  right  hand  side  of  equation  (4.24)  must  toe 
evaluated  at  every  field  point  over  the  entire  wing  surface 
and  summed.  This  sum  will  give  the  street  produced  oy  the 
vorticity  over  the  wing.  Also,  for  each  control  point,  the 
"ight  hand  side  of  equation  (4.39)  must  be  evaluated  at 
each  station  along  the  entire  trailing  edge  of  the  wing  and 
summed.  This  sum  will  give  the  effect  produced  bv  tne 
trailing  vorticity.  When  both  of  these  are  adcec ,  trie 
total  vorticity  effect  at  that  particular  control  point  is 
evaluated.  This  may  be  expressed  as  a  combination  or 


equations  i 4,24.’  and  (4,39)  as  follows, 
total  induced  velocity 


again  two  non-oi  mensi  onal  circulation  -functions.  The 

development  is  similar  to  that  presented  above  tor  the 


circulation  over  the  wing  and  results  in  the  foi lowinc 

A.  'N- 

V^s  Pu  ( y )  -  i-’o0~  p  <Y\>  and  c/s  =  2/AR 

A 

rt(y>  =  (2/AR)P|.0\) 


(4. 2V  ) 
( 4 . 30 ) 


may .  a  m 

<4h/  AR  Ur\7 

'iS\  .  _L  /if?\ 

dij  )  AR  \dr\_  1 


orl  and  dy  = 


<  4.31 


rt ^  =  <xt  -  xP^  +  <n  -  yiP> 


(4.33) 


=  !-[<: 


t  -p 


)  /  r*] 


( 4 . 34) 


Equation  (4.2.)  may  now  be  written  as 


V  = 


G  /d£)  . 

K-  n.r)  v<W  ^ 


(4. 35) 


Now,  transforming  to  angular  coordinates 


pt  (y^>  =  n.  to> 

(lJL'\on .  ms 

W-*J  n  ~  Ue 


Equation  (4.35)  now  becomes 


J_j 

i  j 

r  s  i 

Mrn 

4rr  1 

r  (n-V) 

\6el 

(4. 38 i 
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Equation  <4.24')  is  the  governing  equation  tor  the  induced 
velocity  at  each  control  point  due  to  the  summation  or  the 
streets  c+  every  field  point  on  the  wing  pianform.  Note 
tnat  all  factors  inside  the  integral  are  a  function  of  wing 
pianform  geometry  except  for  the  two  partial  derivatives  of 
the  circulation  function.  Also  note  that  equation  (4.24) 
becomes  indeterminate  when  x  =  x  and  .  This 
condition  occurs  when  the  control  point  and  the  field  point 
coincide  and  is  discussed  later. 

C.  TRAILING  VORTICITY  TRANSFORMATION 

Marty  of  equations  (3.1>  through  (3.  14)  are  identical  at 
the  trailing  edge  and  simply  require  the  subscript  t  to 
be  added  to  the  appropriate  variable.  However,  since  J  is 
constant  at  +1  along  the  trailing  edge  a  few  of  the 
equations  change  slightly  and  are  rewritten  here  for 
c 1 ar i t  y 


i 


■■  4 

.y.ij 


7)  becomes  x ^  = 

X 

r 

+ 

n 

<  4 . 20 ■ 

1 

8>  becomes  Xp  = 

"Lp  +  cp 

<  4 , 2<b  ' 

12)  becomes  dx t 

=  C  h  -  <2/  AR>  2TD(Tdr\_ 

<  4 . 27  ) 

14)  becomes  d;:^. 

=  yUC  dr^ 

<4. 28) 

-  i 

The  circulation  aiong  the  trailing  edge  of  the  wing  is 


a  function  of  the  spanwise 


coordinate  only. 


Cons l der 


from  wmcn  it  foi  lows  that 


tc'-ia:  1  o  n 


i  now  becomes 


O  TT 


.4.2: 


6-VC 


Upon  reversing  the  limits  or  integration  and  re/ersing  the 
algebraic  sign  of  the  integrand,  equation  (4.23)  can  be 
wn  tten  as 


diving  two  equations  which  reduce  to 


Z_  j_  (w 
AR  c  \  <3  \ 


ar  LUn./  c  ui 


(4. 10; 


Now,  the  following  expression  that  occurs  in  equation  (4. 1 


. At) 


+  vv  —  v, 


,w\ 


AL.I 

Zu) 


can  be  expressed  in  terms  of  |  and  as  foli 


4 ,  12) 


Regrouoing  terms  in  equation  (4,12)  and  setting 


F  =  (x 


!p>  - 


dx  dy  =  c  d  diq_ 


i  4 .  13) 


Equation  (4.1)  may  now  be  written  in  the  following  form 


p  4ir\AR, 


,  I 

A{Ffl\)-h  c^‘n^(l£)}c!!drL  <4-15’ 


B.  CHANGE  OF  F  TO  ANGULAR  COORDINATES 

Utilising  equations  (3.16)  and  (3.17)  together  with  a 
non-dimensional  circulation  function  in  the  (^,0) -plane  the 


following  t r ansf or mat l on  can  be  made 


reters  to  values  along  the  trailing  edge  of  the  wing.  It 
:s  desirable  to  transform  these  equations  to  equivalent 
forms  in  the  angular  coordinate  system. 

A.  CHANGE  OF  P  TO  WING  COORDINATES 

Circulation  has  the  dimensions  of  velocity  times 
distance.  Consider  two  separate  circulation  functions,  one 
non-dimensional  with  respect  to  semi -span  and  free  stream 
velocity,  and  one  non-dimensional  with  respect  to  mean 
chord  and  free  stream  velocity.  The  following 

rel at  1 onshi ps  can  then  be  established 

A  /V 

P(x,y)  =  V^c  r<f  ,Y^)S  c/s  =  2/AR  (4.4) 


A  />-/ 

P<x,y>  =  (2/AR)  r<|,T\.)  (4.5) 


Substituting  equations  (3.1)  and  (3.14)  into  ( 4 .  t-i  gives 


rtegrouping  terms 


The  coefficients  of  d|  and  of  drj^must  separately  vanish. 


IV. 


GOVERN  I NG  tiQLiAT  I ONS 


The  vertical  velocity  induced  at  any  particular  canti"oi 
point  consists  of  that  contribution  made  by  the  vortiaty 
distribution  over  the  wing  plan-form  plus  the  contribution 
made  by  the  trailing  vortices  behind  the  wing.  Farrs 
TRe-f,  13  provides  a.  quite  detailed  derivation  of  the  two 
general  equations  +or  induced  velocity.  Both  equations  me¬ 
al  sc  be  found  in  many  other  aeronautical  references  such  as 
Bert  in  and  Smith  CRef.  2  3  and  Kuethe  and  Chow  CRef .  33. 
The  fallowing  farm  of  these  equations  was  utilized  in  this 


and  rt^  is  obtained  by  replacing  with  ^  in  equation 


(4.3)  . 

Equation  (4.1)  is  basically  the  same  as  Park  s  equation 
(3.6)  where  w  '  is  the  velocity  induced  at  a  control  point 
due  to  the  vorticit.y  distribution  over  the  wing  surface. 
Equation  (4.2)  is  the  same  as  Park's  equation  (3.8)  where 
w  ''  is  the  velocity  induced  at  a  control  point  due  to  the 
trailing  vorticity  behind  the  wing.  Subscript  t  here 

IP 


0  =  0  0  =  rr 


Figure  2.  Linear  Coordinates  (1  ,\)  versus 
Angular  Coordinates  (0,5) 
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Trailing  Edge 


ANGULAR  COORD i NATES 


it  is  oesi'anie  to  transform  the  wing  plan-form  from  zhs 
(  $  ,Yi )  —  plane  into  the  ( 0,  6  >  -p  i  ane .  Tms  is  accompiisr 
using  the  relationships 


;nea 


0.5  * 

(  1  -  cos  P  ) 

■m* 

a 

cs*m 

=  0.5 

sin  ^  d 

0 

<  3- 

1 6  > 

cos  & 

;  d  T|_  =  -si  n 

8  dB 

<  3. 

1  7 ) 

at  l  on 

<3. 16) ,  as 

\  van 

es  from  zero 

at  the  1 

ead 

i  no 

1  at 

the  trailing 

edge , 

<0  wi  i 

1  vary 

from  zero 

to 

TT. 

vel  y . 

From  equat 

i  on  ( 3 

.  17)  , 

as  ry^varies  from 

i 

B.TL 

the  left  wing  tip  to  +1  at  the  right  wing  tip,  0  will  varv 
from  TT  to  zero,  respectively.  These  limits  are  important 
in  subsequently  developed  integral  equations.  An  apdeo 
bonus  of  this  transformation  is  that  it  results  in  mesn 
paints  near  the  leading  and  trailing  edges  ana  wing  tips  -:o 
be  closely  spaced,  allowing  finer  resolution  in  tnese 
critical  areas  of  analysis.  Fig.  2  shows  the  relationship 
between  the  linear  coordinates  anP  the  angul  ar 
coordinates  (  0,9). 


£-.  NON-DIMENSIONAL  COORD  I  NATES 

rhe  tirst  change  01-  variables  involves  a  transtormat  1  o 
troa  the  <x,y)— plane  into  the  ( 1  ,f\)  -plane.  In  th 
following  r ei at  1 onshi ps ,  non-subscr i pted  variables  refer  r 
field  points,  while  variables  with  subscript  'p  refer  t 


control  points. 

y  =  >\_  ;  dy  =  df\_  (3.1 

yP  -  nP  >  0 

hl  =  Ain\.  <3.3 

*Lp  =  ?<\?  (3.4 

c  =  (2/AR)  C(1  +T)  -  2  T  <TT\J  (3.5 

Cp  =  (2/AR)  C(1  +T)  -  2'T"^_p  3  (3.6 

x  =  X|_  +  c  \  (3.7 

xp  =  xLp  +  Cp]p  (3.8 

dx l  =  ^•(rdy^  (3.  9 

dc  =  -  (2/AR)  2r<r  dt\_  (3.10 

dx  =  dx^  +  cd 1  +  \dc  (3.  ii 

=  A<Tdirt  +  cd\  +  ^  C  —  ( 2 /  AR )  2 T"<r dr\_  1  (3.12 

/^=Cf\-(2/AR)2'r?3  <3.13 

dx  =  cdj  +jU<Tdr\_  <3.14 

T  —  ( 1  —  taper  ratio)  /  (1  +  taper  ratio)  (3.15 

This  transformation  simplifies  the  algebra  of  th 


problem  significantly.  The  chordwise  dimension  now  vane 
from  0  at  the  leading  edge  to  +1  at  the  trailing  edge 
while  the  spanwise  dimension  varies  from  -1  at  the  1 ef 


wing  tip  to  +1  at  the  right  wing  tip. 


WING  GEOMETRY  AND  COORDINATE  TRANSFORMATION 


The  wing  geometry  considered  in  this  development 
consists  of  a  wing  plan-form  symmetric  about  midspan  witn 
straight  leading  and  trailing  edges  containing  no 
discontinuities,  except  at  midspan  where  allowance  is  made 
for  sweepback.  The  wing  has  no  control  surfaces.  A 
typical  wing  is  pictured  in  Fig.  1.  with  the  geometric 
coordinates  labeled.  The  wing  is  completely  specified  by 
three  factors,  the  aspect  ratio,  taper,  and  the  ieacing 
edge  sweep. 

All  coordinates  on  the  wing  are  normalised  to  a 
semi -span  length  of  one  unit.  Sweepback  at  the  wing  tip  is 
Lamda(A).  EtdO^l  is  the  normalized  spanwise  coordinate. 
Parameter  Tau (T)  defines  wing  taper.  Xi is  the 
normalized  chordwise  distance  from  the  leading  edge. 
Sigma  «T)  is  a  sign  parameter,  +  1  for  the  right  wing  and  -1 
for  the  left  wing.  The  mean  chord  (c)  is  a  function  of 
aspect  ratio(AR),  that  is,  c  =  2/AR.  This  analysis 
involves  two  coordinate  transformations.  The 
transformations  proceed  from  the  (x,y)— plane  to  the 
non-dimensional  <  ^  -p  1  ane ,  then  to  the  anouiar 


bssi c  ana  additional  lift  di stri butions  are  both  known,  the 
c-nrrespondi  ng  induced  drag  can  readily  be  calculateo. 

Inis  analysis  develops  a  metnod  to  calculate  the 
spanwise  and  chordwi se  pressure  distributions  over  the  wing 
which  are  necessary  in  designing  the  wing  structure. 

Although  Parks'  calculations  did  not  achieve  a 
satisfactory  result  ,  he  pointed  out  areas  o-f  concern  for 
possible  errors  and  indicated  a  need  tor  -further 
i n vest i gat i on . 

Where  Parks  used  a  staggered  mesh  between  control 
points  and  -field  points  on  the  wing  to  avoid  the 
indeterminacy  encountered  whenever  a  control  point  and 
field  point  coincide,  the  present  analysis  resolves  the 
indeterminacy  and  works  with  control  points  and  -field 
points  superimposed  upon  each  other.  Also,  the  series 
solution  employed  by  Parks  to  represent  the  circulation 
distribution  has  been  discarded.  Instead,  this 
distribution  is  now  represented  more  simply  and  directly  by 
the  values  o-f  the  circulation  at  the  field  points 
themsel ves . 


I N T RuDUC T I UN 


1 1 


This  analysis  is  a  continuation  of  a  previous  siu.ov 
conductec  by  John  L.  Parks,  Lieutenant  Commander ,  ur 1  tea 
States  Naval  Reserve  [Ret.  1]. 

The  concept  or  a  continuous  vortex  sheer  of  van  sole 
strength  over  the  wing,  trailing  to  infinity  aft  or  tne 
wing  is  utilized  to  determine  the  value  of  the  circulation 
function  at  a  finite  number  of  control  points  on  tne  wing. 
The  vortex  sheet  strength  is  restricted  by  tne  r-utta 
condition  at  the  trailing  edge  of  the  wing  ana  by  the 
requirement  of  no  flow  through  the  wing  at  the  specif iea 
control  points.  The  present  analysis  is  restricted  to 
steady,  inviscid,  incompressible  flow  about  thin  wings  witn 
straight  or  swept  leading  and  trailing  edges. 

Wing  aerodynamics  involves  two  fundamental  aspects, 
first,  a  desiren  basic  lift  distribution  may  oe  specif  iso 
over  the  wing  and  the  correspond! ng  wing  caraoer  must  be 
determined.  Second,  a  flat  plate  wing  may  be  specifier 
with  the  resultant  additional  ii+t  distribution  to  ds 
calculated.  Once  the  additional  lift  distrioution  is 
known ,  other  pertinent  aerodynamic  parameters  may  oe 
Determined;  notably,  the  slope  of  the  wing  lift  curve  anq 


the  location  of  the  aerodynamic  center 


Also,  when  the 
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Equations  (4.24)  and  *4. 39)  may  be  evaluated  by  several 
ait+erent  methods.  The  relationship  between  the  geometric 
placement  of  the  control  points  and  tieid  points  is  parnaos 
one  ot  the  most  important  factors.  wing  symmetry  scout 
mi  dspan  permits  the  analysis  to  utilize  control  points  eve'" 
only  one  semi— span  of  the  wing.  The  e+f ect  ot  field  points 
and  trailing  vorticity  must  include  the  entire  wing 
pi anfocm:  however ,  symmetry  about  miespan  again  permits 
utilization  o+  numerical  values  over  only  one  semi -span. 
The  right  semi -span  is  chosen  for  this  analysis. 

As  mentioned  before,  equation  (4.241  becomes 
indeterminate  if  the  control  point  mesh  and  the  field  point 
mesh  coincide.  Equation  (4.39)  also  becomes  indeterminate 
when  the • spanwi se  coordinate  of  the  trailing  vorticity  ana 
the  concerned  control  point  coincioe.  Note  that  S  in 
equation  (4.34)  approaches  zero  as  the  spanwi se  coordinate 
of  the  control  point  and  the  trailing  vorticity  converge. 
One  method  to  avoid  these  indeterminate  points  is  to 
stagger'  the  different  meshes  so  that  they  never  coincide. 
Another  method  is  to  determine  an  analytical  limit  or 
principal  value  of  the  integral  as  the  indeterminate  point 
is  approached.  The  second  method  is  chosen  for  tn  = 
analysis  and  is  presented  in  the  following  paragrapns. 
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uihers  a  concroi  point  and  a  field  point  coincide,  me 
numerator  and  the  denominator  or  equation  (4.24;  are  bctn 
cero  resulting  in  an  indeterminacy.  Since  an  i ndetermi nacv 
can  only  occur  at  a  control  point,  all  occurrences  are 
con-fined  to  the  right  semi— span.  All  geometric  variables 
ana  sigma,  the  sign  parameter ,  have  a  positive  value  on 
this  section  of  the  wing.  Fig.  3  depicts  an  exploded  view 
or  a  mesh  point  singularity.  The  discrete  value  or  the 
circulation  at  the  singular  point  must  be  determined. 
Recal  1  that  the  circulation  -function  itsel-f  is  continuous 
over  the  wing;  therefore,  it  is  continuous  m  the  proximity 
of  the  singular  point.  Variables  '  Lt '  and  'v  are 

introduced  to  represent  small  deviations  in  <f>  ana  0.  in  the 
vicinity  of  the  singularity, 
contribution  of  the  singular 
velocity  at  the  singular  point, 
i s  estabi i shed 

+<f  +( 

AR  '  = 

where  symbol  W  is  an  abbreviation  that  stands  for  the 
integrand  of  equation  1 4 . 24  > .  Now,  let 

u  =  (  i  5 . 2 i 


If  AR ‘  represents  the 
element  to  the  induced 
the  following  reiaticsMp 


^  / 


Figure  J.  Field  Element  at  Singular  Point 


equat l on 


where  F'  represents  the  expression  in  brackets  in 
(4.24)  . 

An  analysis  of  equation  (4.24)  in  the  irrimediate 
vicinity  of  the  singularity  results  in  the  following 

+  (2/AR )  l  ( 1  +T>  -  1  (5.5> 

xp  =  Aqp  +  <2/AR)C(l  +  T>  -  2^  1  lp  (5.6) 

<x-xp)  =  A(r\-qp)  +  (2/AR)  C  <1+T)  (?-ip)-2T(r\?-«lp?p)  1  (5.7; 

rx  =  cos  9  =  cos(0p+v)  =  cos  —  sin  5p  v  (5.8) 

%  =  cos  £p  (5*9) 

<l\  -  qp)  =  -  sin  0p  V  (5.  10) 

\  =  0.5(1  —  cos  (j))  -  0.5C1  —  cos  (^p+u)  1 

=0.5(1  -  cos  ^  +  sin^  u)  (5.11) 

]p  =  0.5(1  -  cos  0p)  (5.12) 

(  |  -  Ip)  =  0.5  sin  ^p  u  (5.13) 

“  rlp?p)  =  (rlp-sinep  v)  (]p+0.5sin^p  u)-nplp  (5.14) 

=  qp"ip  +  nP  0.5sin^p  u  -  Sin6p  v  ?p 

-  sin  6 p  v  0.5  sin  0p  u  -  r\p]p  (5.15) 

Simplifying  and  neglecting  higher  order  terms  in  ,  v~ , 


and  uv 


Oxl-ir^plp)  =  C  (rjp/2 )  sin  <f>  p3u  +  (-sin  0p  ?p)v  +  ..-(5.16) 

Substituting  equations  (5.10),  (5.13),  and  (5.16)  into 

(5.7)  , 

(;•;  -  x  )  *  Au  +  Bv  (5.17) 

where  A  and  8  are  known  -functions  o-f  0p,  0p  ,  *\p ,  J  p  ,  T? 
and  AR.  The  denominator  of  the  integrand  o-f  equation 
(4.24),  in  the  immediate  vicinity  o-f  the  singular  point  may 
now  be  expressed  as 

r -  C  (Au  +  Bv)^  +  (-sin  0p  v)2l^"w'  (5.18) 

Since  6n  is  -fixed,  let  -sin  0p  =  C,  then, 

r3  =  C (Au  +  Bv)2  +  C2v2l 1 • 5  (5.19) 

Assume,  by  analogy  with  the  above  development,  that  the 
numerator  o-f  equation  (4.24)  may  also  be  expressed  as  a. 
■function  o-f  u  and  v,  namely 

F’  =  Du  +  Ev  (5.20) 

Now,  switching  to  polar  coordinates  with  '  the  polar 
radius  and  '  t*  '  the  angle,  let  u  =  p  cos  at  and  v  =  sin  cL  . 
The  effect  o-f  the  singular  field  element  must  be 
determined.  An  integration  over  the  entire  element 


produces  the  required  result. 


The  integration  covers 


total  angle  of  21 T  and  a  polar  radius  varying  -from  zero  at 
the  center  of  the  element  to  its  outer  bounds.  Equations 
(5.19)  and  (5,20)  expressed  in  polar  coordinates  become 

r°  =  C  (Acos  <*.  -Bsi  n  ol)  si  n^  <*■  3  *  ■  w'  (5.21) 

=  J  <K)  p  *'  (5,22/ 

where  J  (®0  is  an  abbreviated  -form  of  the  term  in  brackets. 


F'  =  CD  cos  «.  +  E  sin  oe.  ] 


=  H(aO 


(5. 23) 


(5.24) 


where  H  (<*•)  is  also  an  abbreviated  form. 

Note  the  polar  symmetry  of  each  of  the  above  functions, 
that  is 


J  (  o<.  +■  VC  )  =  J  (oL) 


(5.25) 


H(Qt  +  tT)  =  -H(oC)  (5.26) 

Therefore,  if  &R  '  is  the  effect  of  t  ne  1 ndet ermi nacy 


AR'  = 


f  fjlfc')*  ( 

J  J  GWe3 


dc<  d 


c  !  = 


(5.2?) 


The  value  of  this  integral  is  zero  by  reason  of  odd 
polar  symmetry.  This  proves  that  the  effect  of  the 
singular  field  element  on  the  control  point  is  zero  and  mav 
be  neglected. 


B.  INDETERMINATE  POINT  ALONG  TRAILING  EDGE 

As  approaches  *r^_.  in  equation  <4.39)  ,  G  a.pproache 
zero;  therefore,  the  integrand  is  i ndeter mi nate .  Note  tha 
(x^.  -  Kp)  does  not  approach  zero,  since  Xp  is  th 
x  —  coordi  nate  of  the  center  of  the  concerned  control  poin 
and  Xi  is  the  x-coordinate  at  the  trailing  edge  of  th 
wing.  In  the  i mmedi ate  vicinity  of  the  singularity 
( —  Y|_p )  becomes  very  very  small,  but  remains  finite 
therefore,  it  can  be  represented  as  a  small  value  times 
constant.  Let  this  constant  be  ( -  Xp)  ;  therefore,  a 
( -  T^p )  approaches  zero,  the  following  relationships  ma 
be  established 


( 5 .  28 


<  5.  29 


r«\.  -  )  =  e*<xt  -  Xp)  ,  where  e  <<  1 


V  =  <>:t  -  V"  +  <«\  - w 


=  (xt  -  xp)‘(l  +  e-) 


rt  =  ( x t  -  xp)  <  1  +  e^  )  -  • 

=  ( x  t  -  xp)  (1  +  e2/2  +  ...H.O.T.) 


The  integrand  of  equation  (4.39)  may  now  be  expressed  as 


IV  ’IP) 


e(V*p)  (.  (VW  + 


which  can  be  reduced  to 


(h  ~  x?)  \  I  +  e% 


Thus,  in  the  limit,  as  ^  approaches  "*\p  ’  e  approaches 
zero,  and  the  value  of  equation  < 5 . 32 )  becomes  zero.  This 
proves  that  the  trailing  vorticitv  contribution  directly 
downstream  of  a  control  point  makes  no  contribution  to  the 
velocity  induced  at  the  control  point. 


VI 
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ANALYSIS 


Consider  a  given  wing  plan-form  for  which  the  three 
parameters  aspect  ratio,  taper ,  and  leading  edge  sweep  are 
known  values.  The  correspond! ng  lift  distribution  and 
therefore  the  circulation  distribution  over  this  wing  are 
quantities  to  be  determined.  Ail  quantities  in  the 
integrands  of  equations  <4.24)  and  <4.39)  are  known  except 
for  the  partial  derivatives  of  the  circulation  function. 

A  mesh  site  adequate  to  represent  the  desired  lift 
distribution,  while  not  requiring  an  excessive  amount  of 
numerical  caicul ations,  is  desired.  A  mesh  consisting  of 
eight  chordwise  stations  and  sixteen  stations  across  the 
entire  wing  span  is  arbitrarily  selected  and  should  produce 
the  desired  accuracy.  As  previously  mentioned,  wing 
symmetry  about  midspan  allows  all  computations  to.  be  made 
with  the  numerical  values  on  the  right  semi-span. 

A.  WING  GRID  POINTS 

This  analysis  involves  mapping  points  on  the  wing  from 
the  physical  <x,y)— plane  to  corresponding  points  in  the 
‘~4>,6'>  —plane  and  vice  versa.  Fig.  4  depicts  a  uniform 
rectangular  grid  of  points  in  the  lane,  constructed 
as  foil ows 

0.  =  40fi-O.5>;  40  =  26  =17/ 8;  i  =  1,2,...,  8  <6.  1 

34 


9  ,  =  A6  (  J  — 0 . 5  )  ;  A0  =  2  d  =  TT  /  1 6 ;  j  =  1,2,.  ..,8  (6.2) 

Also,  along  the  leading  and  trailing  edges  of  the  wing 

0L  (  j  )  =  0;  <£>t  <  j  )  =  IT  (6.3) 

while  along  the  wing  tip  and  at  midspan 

0t<i>  =  0;  0m(i)  =  TT/2  (6.4) 

The  above  equations  will  produce  64  grid  points  on  the 
right  semi-span  with  each  point  centered  in  its  respective 
rectangular  mesh  block.  The  'j'  subscript  varies  spanwise 
-from  right  wing  tip  to  midspan  while  the  'i  '  subscript 
varies  chordwise  from  leading  edge  to  trailing  edge.  The 
unknown  values  of  the  circulation  function  at  each  of  these 
grid  points  must  be  determined. 

For  every  field  element  on  the  right  wing  there  exists 
a.  c.orrespondi  ng  mirror  image  field  element  on  me  i  ctt 
wing.  By  utilising  the  symmetry  between  each  fieid  element 
and  its  mirror  image,  the  range  of  the  'j'  subscript  can  d? 
restricted  to  8  vice  the  16  it  would  otherwise  be.  Fig.  55 
depicts  the  planform  distribution  of  the  circulation 
function  with  appropriate  subscripts  for  matri:;  and 
corresponding  vector  represent at i on ,  All  quantities  with 
numerical  values  at  respective  grid  points  may  pe 
designated  similarly.  Should  it  become  necessary  to 
transform  from  matri;-;  notation  to  vector  notation  or 
vice  versa  the  following  rel at i onsh i ps  are  used,  with  k 
denoting  the  respective  vector  subscript. 

k  =  SCi  -  1>  +  j  where  Pv v k )  =  < j ,  i  )  (6.5) 
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Figure  5.  Gamma  Matrix  and  Vector  Grids 
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1 


l:;  /  H  +  1 


■  k/y  by  integer  division) 


i  a .  6 ) 

j  =  k  -  8  (  i  -  1  >  where  (  j  ,  l  )  =  H,  ( k  >  <6.  7 ) 

B.  VARIABLE  NOTATION 

The  subscript  notation  introduced  in  the  aoove 
paragraphs  is  used  extensively  in  the  subsequent  analysis. 
The  -following  notation  -for  various  functions  and 
derivatives  is  also  used.  Note  the  different  expressions 
that  correspond  to  the  same  numerical  value. 

The  first  derivatives  of  the  circulation  function 


rhe  first  derivatives  at  the  leading  and  trailing  edge 
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(J)i  1  =  <  1  /bj>  >Li-l  /  2 )  n  i  j  ,  i  -  1  >  i-  <  1  /  2  ?  P  \  j  ,  i  f  i  i  j 

r^(j,3;  =  il/b4»  C  <-3/13>rc  j  ,  7)  +  <  3/ 13)  P\  j ,8?  1 


A1  i  values  on  the  Is+t  hand  sides  of  the  abo'- 


equations  are  determined  by  equation; 


:■)  and 


i he  solution  tor  the  circulation  values  requires  inversion 
o+*  a  3  by  S  matri:; . 

Wow,  equations  (7.17),  (7,18),  and  (7.15?  are  utilized 
to  calculate  the  spsnwise  derivatives  of  the  circulation 
function.  All  variables  or,  the  right  hand  sides  of 
equations  (7.24)  and  (7.25)  are  now  known.  The  wine  slcpe 
function  may  be  written  as 


‘V  -  I 

k 


where  S  is  a  64  by  64  matri;-:  of  influence  coef  f  ici  ent- 
whose  e-:act  definition  is  implied  by  equations  (7.24?  an; 
(7,25?,  A  more  detailed  analysis  of  matri'.;  B  is  given  : 
Ref .  4 , 


A.  ADD  1 T 1 CNAL  LIFT 

In  the  additional  lift  problem,  wp  =  -1, 

pressure  distribution  over  the  wing  is  desired, 
following  previously  developed  equations  are  reor i rec 

<  j  ,  i  >  =  a  >  c  r<  j ,  i )  +  <  i  /3)  r<  j  ,2?  j 


;np 


,7.9) 


i 


a^i  ■■ 


Equations  (7.2) ,  <7.24)  and  (7. 25) ,  when  evaluated  at 

each  control  point  over  the  *~sp?n ,  produce  si^tv-^OLir 

equations  in  the  sixty— tour  unknown  values  of  the 
circulati on  t uncti on .  The  sol ut i on  to  this  system  ot 
equations  produces  the  desired  values  ot  the  circulation 
function  at  the  mesh  points  over  the  right  semi— span. 

E.  BASIC  LIFT 

In  the  basic  litt  problem  a  pressure  distribution  is 
soecified  over  the  wing,  leaving  the  corresponding  wing 
slope  function  to  be  determined.  In  this  case,  equation 
<6.225  ,  rearranged,  gives  the  chordwise  derivatives  o+  the 
circulation  function  at  each  mesh  point  on  tne  wing,  as 
f  o  1  3  ows 

^  C  (-j)  Sua  P(^x)  <7.26) 

At  the  leading  edge,  equation  <6.28)  applies 

r*L<j>  =  ^  pL(j)  <7.27) 

Now,  equations  <7. 10) ,  <7. 14) ,  and  <7. 16)  are  expressions 

for  in  terms  of  the  unknown  circulation  values.  These 

equations  are  repeated  below  with  known  values  collected  on 
the  left  side  of  each  equation. 

fyu,i>  +  <  a  <p/  4 )  r^L  <  j  >  =  <  i  /a0)  6  r<  j ,  i  >  <  7 .  i «.’  > 


JjjL 


.  F  t  ti_D  tLtiitNT  ANl>  i  MeGt  ELEMEN i 

Pec-si  1  rhst  tre  oni  v  d;  Aorencss  in  the  values  ove 


i  g t- f  and  i  e  +  t  semi  -soan 


sro  the  algebraic  signs  of 


nd  CT.  Therefore,  w  ’  may  be  thought  of  as  having 
acts,  one  pact  consisting  ot  the  contribution  th 
articular  tield  coin!  on  the  riaht  semi— scan  males 


specif  isa  ccnccoi  point  ,  ana  anorr.er  pare  consisting  nt  rn 
contribution  that  the  image  field  point  on  one  let 
semi— span  raves  on  tr.e  same  control  point.  The  5  am 
concept  applies  to  ,  the  trailing  vortici tv  inteaca 
along  the  trailing  edge.  Eauations  <4, 2d)  and  { d. 3^ >  ar 
mod  1 f i ed  sliohtly  so  tnat  the  eff ec t  of  a  tie! a  point  o 
the  right  semi— span  and  the  effect  of  its  1  mage  on  the  1?+ 
semi  -span  may  be  calculated  si  mu  1 taneousl y.  us: no  symmecc 
relations.  Let  the  subscript  i  represent  the  riqn 


semi -spar.  ana  1  et  supeccipt 


cepresent  t.ne  i  =t 


~emi  -span  .  -'la  1  >g  rr-c  s>  I  tutions  f  ijr  the  =jiiiT,crr’ 


--me1-,  t  c. 


a  t  i  :•  n  a  ;  4 . 2  4  >  r.  d  •  4  . 


.n  3  /  (p  r  tp\. 


cor  tii  cions  may  be  used  in  this  case.  tor  the  fignt  wing 
tip  !T03n  points 

^  *  1  f  i  )  =  ( 1  /AS >  [  P(  1  ,  i  >  +  <  1/3)  P(2,  i  >  3  '  7 .  i  7  1 

Central  difference  formulas  are  also  the  same 

re<j,i>  =  <  1/A0)  L  (-1/2)  T(  j-l  ,1  )  +  ( 1/2)P(  j  +  l  ,1  )  ]  (7.  18) 

At  midspan,  allowance  is  made  for  a  discontinuity  in 
the  slope  of  the  circulation  function  as  mentioned 
previously.  Again,  a  polynomial  suitable  to  satisfy  the 
boundary  conditions  is  assumed;  therefore,  let 

(T-Toi  =  a,  <  0  -  9o>  +  < i / 2 > a.-> (6  -  0 o 1  ^  (7,i?) 

uj  i.  O  i— > 

and 

^  =  At  +  A2f0-0g)  (7.20) 

Satisfying  equation  (7.19)  at  j  =  6  and  j  =  7,  and  solving 
for  A,  and  A-,,  gives  the  needed  coefficients  in  equations 
(7.19)  and  (7.20).  Now,  applying  equation  <7. 19)  at 
midspan  and  equation  (7.20)  at  midspan  and  at  j  =  8  gives 

rm  (  i  )  =  ( 3 / S >  P < 6  ,  i  >  -  ( 5 / 4 >  P ( 7 , 1  )  +  ( 1 5 / 8 ) P < 8  ,  i  )  (7.21) 

rdm  (  i  >  =  ( 1 /A9)  C  T(6,  i  )  -  3  r  (  7  ,  i  )  +  2p8,i)]  (7.22) 

Pq(8,  i  )  =  (  1  /A©)  C  (  1/2)  P(6,  l  )  -2  P< 7 ,  l  )  +  (3/2)  He,  i  )  J  (7.23) 


V(L  =  -  >  1  2  >  A  ^  <  'TT  -  <f>  ) 


Di  f  f  erent  1  at  1  ng  again, 


=  +A-T  (It  -  ^ 


Applying  equation  (7.11)  at  two  trailing  grid  points, 
solving  tor  A-r ,  and  substituting  into  <7.  12)  ,  a  di  ffei-ence 


expression  for  1^  in  te 


rd  in  terms  of  !*(  j  ,  7)  and  P(  j  ,8)  is 


Obtai  nprj 


as  follows 


<  1  /A0)  C  (-3/ 13)  P(  j  ,  7) 


(3/13)  r ( j ,8) 


>  J  (7.1. 


A1  so 


rt(j>  =  (-1/26)  P  ( j  ,7)  +  (27/26)  Pj  ,8) 


For  mesh  points  2  through  7  in  the  chordwise  direction, 
central  differences  are  used  which  results  in  the  following 

<  j  T  i  )  =  (  i/Atf>>  C  (-1/2)  P<  j  ,  i-1 )  +  (l/2>P(j,i+l)  3  (7.16) 


The  difference  equations  for  the  trailing  edge  an.n 
central  points  are  al  so  valid  f  or  the  basic  lift  prctsi 


C,  SF'ANWISE  DERIVATIVES  OF  CIRCULATION 

Along  the  right  wing  tip,  the  circulation  is  zero: 
however,  the  first  derivative  is  finite.  these  conditions 
are  precisely  analogous  to  the  boundary  conditions  along 
the  leading  edge  for  the  additional  lift.  Hence,  the  same 
difference  ecuation  developed  above  for  these  boundary 


Equations  (7.6i  and  (7.7)  satisfy  the  boundary  conditions 


■for  the  additional  lift.  These  are  that  F  vanish  and  F^  is 
finite  at  the  ieacing  edge  where  -  0.  Applying  equation 

(7.6)  at  two  grid  points  adjacent  to  the  leading  edge, 
solving  for  A,  and  A-,,  and  substituting  into  equation 

( 7.7)  ,  a  difference  expression  for  F^  in  terms  of  T j  and  f7 
is  obtained  as  follows 


=  <  1  /A 4> >  C  F,  +  ( 1 73 )  V21 


Therefore,  for 
the  following 


any  mesh  point  adjacent  to  the  leading  edge, 
subscript  notation  is  used  for  Y ^ 

( 1/A0)  C  P<  j  ,  1 )  +  <  l/3>r<  j  ,2)  3  (7.9) 


A  similar  analysis  applied  to  the  basic  lift  results  in  the 
foil owi ng 


=  (1/A^)  C-(A0/2)2P^l<  j)  +  6r<j,l):  (7.10) 

The  procedure  for  determining  the  difference  formula 
applicable  near  the  trailing  edge  requires  defining  a 
function  for  the  circulation  that  is  numerically  equal  to 
|”y  at  the  trailing  edge  with  first  and  second  derivatives 
equal  to  sera.  This  requires  beginning  with  a  cubic 
equation;  therefore,  let 


P-rt 


(  1  76)  At  ('TT  -  4>  )  3 


1  1  ) 


Di f f erent l at l ng  once, 


tx  or essi  on 


s  f nf  the  partial  derivatives  of  the  circular 


function  m  ter^s  of  the  unknown  circulation  values  ar  = 
raeded  and  are  derived  below. 

When  solving  the  basic  lift  problem  with  a  specif i ec 
pressure  distribution,  the  unknown  is  the  wing  si  one 
function.  Wp.  This  solution  involves  multiplication  Py  a 
64  Py  64  matrix.  Put  inversion  of  the  matrix  is  nor 
requi r sa . 


B.  CHORDWISE  DERIVATIVES  OF  CIRCULATION 

In  devel op ing  finite  diff er ence  squat i ons ,  a  poi ynomi ai 
is  passed  through  a  number  of  points  in  the  particular  area 
of  interest.  The  polynomial  used  in  each  case  must  be  of 
sufficient  order  to  ensure  that  the  required  derivatives 
exist,  and  it  must  satisfy  the  boundary  conditions.  The 
finite  difference  equations  utilized  for  the  central 
difference  formulas  require  a  point  on  either  side  of  the 
point  in  question.  Hence,  the  central  difference  equations 
are  valid  for  all  mesh  points  except  those  that  lie 
adjacent  to  the  leading  and  trailing  edges.  For  the 
leading  edge,  let 


r  =  Aj<£  +  <1/2)A2<2>2  +  0C63j 


Diff  erent i at i ng , 


=  A.  +  A  ^6  +  OCfe: 


(7.7) 


sum  ot  two  components  associated  with  the  basic  i 1 f t  ano 
the  additional  lift,  respectively.  Thus 


wo  “  woB 


<oC-  «*B)wpA 


where  wpB  reflec't3  the  slope  of  the  wing  camber  surface 
according  to  the  relation 


wpB  “  01 B 


(7.6. 


whi  1  e 


wpA  “  1 


reflects  the  uniform  change  in  wing  slope  pe'  radian  r?  + 
additional  angle  of  attack. 

Note  that  the  right  hand  side  of  equation  (7.2; 
consists  of  the  sum  of  equations  <4.24)  and  <4.39).  When 
solving  for  the  additional  lift,  the  only  unknowns  in 
equations  <4.24)  and  <4.39)  are  the  derivatives  of  the 
circulation  function.  These  derivatives  are  expressed  in 
terms  of  the  initially  unknown  values  of  the  circulation. 
With  sixty-four  unknown  values  of  circulation,  sixty-four 
equations  a^e  required  for  a  solution.  Each  control  point 
generates  an  equation  involving  terms  of  the  unknown 
circulation  distribution.  This  produces  sixty-four 
equations  in  sixty— four  unknowns  which  must  be  solved 
simultaneously  to  determine  the  circulation  values.  This 
solution  involves  the  inversion  of  a  64  bv  64  matrix. 


i 
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VII.  GENERAL  i-RQBLEM  ANALYSIS 

A.  PROBLEM  SPECIFICATION 

As  mentioned  previously,  tne  calculation  sequence 
involve  determining  the  wing  slope  function  from 
specified  basic  lift  distribution,  or  determining  the  1 
distribution  from  a  specified  wing  slope  function, 
interest  here  are  both  the  basic  and  additional  1 
distributions. 

The  values  on  the  right  wing  span  are  utilized  in 
calculations.  The  right  wing  span  is  divided  ir'c  ei 
chorcwise  and  eight  spanwi se  stations,  resulting 
sixty-four  control  points  over  the  semi— span.  To  sati 

the  condition  of  no  flow  through  the  wing  surface, 
induced  velocity  must  balance  the  normal  component  of 
free  stream  velocity,  or 


l  ) 

voO 


.n  ol 


=  lz/*> 


Setting  the  free  stream  velocity  to  one 
approximating  sin  <X  with  for  small  angles 

g  i  ves 


un  1 1 
of  atta 


<oz/^x>  —  V  =  w_  =  w  '  +  w_"  <  7 

P  P  P 

By  analogy  with  equation  (6.32),  the  overall  wing  si 


as 


function  a 


each  control  point  may  be 


represented 


Again  note  that  the  specific  additional  lift 
coefficient  CL^  is  identical  with  the  lift  curve  slope 
(dCj_/d«0  and  is  numerically  equal  to  the  value  that  the  net 
additional  lift  would  reach  if  linearly  ex tr apol ated  to  an 
additional  angle  of  attack  of  one  radian.  Once  again, 
however,  the  additional  angle  of  attack  ( °<  -  °<g)  is 
restricted  to  small  values. 

Basic  and  additional  lift  are  also  discussed  in  Ref.  I-. 
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function  with  respect  to  yZ>  at  zero  also. 


T  o  summar i z 


eauar  x  on  f  or  m 

AcPit  =  =  °  < 


hex  t , 

consi 

der  the 

condition  that 

the 

P  re 

coef  f i ci ent 

1  s 

infinite 

at  the  leading 

eaqe 

for 

add i t i onal 

i  i  + 1 . 

At  the  I 

eading  edge ,  0  =  0 

and  si 

n  <fi 

therefore , 

the  denominator 

of  equation  (6. 

22)  l 

5 

zero.  Thus 

is  finite 

at  the  leading 

edge 

*  or 

additional  lift.  Summarised  in  equation  -form 


Finally,  f or  the  basic  li-ft  the  pressure  coefticiei 
finite  at  the  leading  edge.  A  similar  analysis  to  the 
chs  trailing  edge  using  L'Hopitai  s  rule  results  in 
second  derivative  of  the  circulation  function  with  re1 
to  0  being  a  finite  value.  Suininarired  in  equation  for; 


Or,  in  the  notation  introduced  in  section  B 


P  (  i ,  i  >  = 


(  8) 

r 

1  1 

1  C.C-0  -oin.  j 

(  6  .  2  2 ) 


D.  LEADING  AND  TRAILING  EDGE  REQUIREMENTS 

The  boundary  conditions  pertaining  to  equation  to. 22) 
for  the  basic  or  ideal  lift,  and  the  additional  lift, 
determine  the  nature  of  the  circulation  function  and  its 
derivatives.  The  Kutta  condition  fixes  the  pressure 
coefficient  at  zero  at  the  trailing  edge.  For  the 
additional  lift,  the  pressure  coefficient  is  infinite  as 
the  leading  edge.  However,  for  the  basic  lift,  the 
pressure  coefficient  at  the  leading  edge  is  finite.  These 
specifications  together  with  equation  (6.21)  or  (6,22; 
determi  ne  _the  nature  of  the  circulation  function. 

First,  consider  the  Kutta  condition  which  requires  that 
the  pressure  coefficient  be  zero  at  the  trailing  edge.  The 
denominator  of  equation  (6.22)  is  also  zero  at  she  trailing 
edge,  since  $)  =  Tf  and  consequently,  sin  =  0.  This  means 
that  ^  must  not  be  finite  or  the  quotient  becomes 
infinite.  Therefore,  f""^  at  the  trailing  edge  must  be  zero, 
resulting  in  an  i ndetermi nacy .  Using  L'Hopital  's  rule  to 
determine  the  limit,  which  must  be  zero,  the  denominatior 
becomes  cos  Tf  =  —  1  and  the  numerator  must  be  zero.  This 
fixes  the  value  of  the  second  derivative  of  the  circulation 
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F*<j,D  =  c  <-i/2)r<j,i-u  +  <1/2)P( j  ,1  M  >  J 

r^<j,8>  =  (  1  '&$>'■  C  <-3  /  1 3 )  P(  j  ,75  +  ( 3  /  13)  P  •  j  ,  8 '  j 
rg(l,i>  =  i  1/A©)  L  Tci  ,  1  >  +  (l/3>r(2,i  )  3 

=  <i/Ae>  c  <-i/2>n  j-i,i)  +  <i/2)r<  j+i  ,1  >  i 
r0<8,x  >  =  <!/&&•>  C  (l/2)F(S,i  >  —2  Pi  7 ,  i  >  +  (3/2)  IF  3,  i  >  ] 


rt<j>  =  (-1/26)  T<j  ,7)  +  (27/26) F  ( j  , 


8) 


The  wing  slope  function  in  this  case  may  be  written  as 


wD<V  =  Z  A<kpfk>r<k)  <7. 

k 

where  A  is  a  64  by  64  matrix  of  influence  coefficie 
whose  exact  definition  is  again  implied  by  equati ons  (7. 
and  (7,25).  Matrix  A  must  be  .  tverted  in  this  case 
determine  the  circulation  function. 

Once  the  circulation  function  is  known,  its  derivat 
with  respect  to  &  is  found  by  substituting  the  known  vai 
into  equations  (7.9),  (7.14),  and  (7.16).  Equation  fg. 

is  then  utilized  to  determine  the  pressure  distributi 
Other  aerodynamic  parameters  may  then  be  readily  caicula 
for  the  wing  in  question. 

The  matrices  and  calculation  procedures  presented 
the  above  paragraphs  are  described  in  detai i  m 
concurrent  analysis  conducted  by  E.  M.  Barber  [Ret.  41. 
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VIII.  SUhMANV 

The  solution  of  the  system  of  equations  developed 
this  analysis  is  best  suited  for  processing  by  mat 
methods  on  a  digital  computer.  A  brief  summary  of 
procedural  steps  required  in  the  calculations  for 
additional  lift  and  the  basi  c  lift  are  given  in 
to! ! owing  paragraphs. 

A.  ADDITIONAL  LIFT 

Given  wn  =  -1: 

Fir st ,  calculate  matrix  A  utilizing  equations  >7. 

( 7 .  1 4  >  ,  (7.15),  <7.16),  (7.17),  (7.18),  and  (7.23)  toget 

with  the  planform  parameters  in  equations  (7.2 4) 
(7.25). 

Next,  invert  matrix  A  and  multiply  by  vector  wp 
determine  the  circulation  values. 

Now,  determine  the  chordwi se  derivative  of 
circulation  function  utilizing  equations  (7.9),  (7.14), 

(7.16). 

Finally,  use  equation  (6.22)  to  calculate  the  press 


i  n 

r  i  'i 

t  be* 
the 
the 


9)  * 
her 
and 

to 

the 

co  r !  b 

u  e 


d i st  r i but l on . 


B. 


BASIC  LIFT 


Gi ven  P ( j  ,  i >  : 

First,  calculate  the  chordwise  derivative  of  the 
circulation  function  using  equation  (6.22) ,  and  the  second 
derivative  factor  using  equation  (6.28), 

Next,  compute  the  circulation  function  using  equations 
(7.9),  (7.14),  and  (7.16). 

Now,  compute  matrix  B. 

Finally,  multiply  matrix  B  by  the  circulation  vector  to 
obtain  Wp. 


T  ’  ,  w.m  .  "l  "'I  ■■l" 


IX.  CONCLUSIONS 

Although  no  numerical  results  are  yet  available,  there 
are  ample  grounds  tor  confidence  in  the  validity  o+  the 
calculation  method  presented  in  this  thesis.  All 
derivations  are  based  on  sound  mathematical  procedures  and 
assumptions.  This  analysis,  along  with  the  parallel  effort 
presented  in  Ret,  4,  has  satisfied  the  initial  purpose  ot 
the  overall  research  effort,  which  was  to  develop  an 
adequate  mathematical  model  for  calculating  the  aerodynamic 
performance  of  a  wing  of  arbitrary  planfonn  in  inviscio. 


i ncompressi bl e  flow 
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